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Abstract. We study the influence of an aperiodic extended surface perturbation on the surface critical behaviour of 
the two-dimensional Ising model in the extreme anisotropic limit. The perturbation decays as a power k of the distance 
I from the free surface with an oscillating amplitude Ai — ( — where /; = 0, 1 follows some aperiodic sequence with 
an asymptotic density equal to 1/2 so that the mean ampltitude vanishes. The relevance of the perturbation is discussed 
by combining scaling arguments of Cordery and Burkhardt for the Hilhorst-van Leeuwen model and Luck for aperiodic 
perturbations. The relevance-irrelevance criterion involves the decay exponent k, the wandering exponent uj which 
governs the fluctuation of the sequence and the bulk correlation length exponent v. Analytical results are obtained 
for the surface magnetization which displays a rich variety of critical behaviours in the (k, uj) -plane. The results are 
checked through a numerical finite-size-scaling study. They show that second-order effects must be taken into account 
in the discussion of the relevance-irrelevance criterion. The scaling behaviours of the first gap and the surface energy 
are also discussed. 
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1 Introduction 



At the free surface of a homogeneous d-dimensional sytem 
with short-range interactions displaying a bulk second-order 
phase transition, the scaling dimension x^ of the surface energy 
density is equal to d [|lj] . As a consequence, a weak short-range 
surface perturbation AK S of the reduced interaction K = (3J= 
J/k-oT cannot change the surface critical behaviour, since its 
scaling dimension dak = d — 1 — x s c = —1 and the perturbation 
is irrelevant. Such a pertubative argument does not exclude the 
occurence in d > 2 dimensions of special, extraordinary, and 
surface transitions for strong enough enhancement of the sur- 
face couplings [||]. 

The situation is somewhat different for the Hilhorst-van 
Leeuwen (HvL) model Q], in which the surface perturbation 
extends into the bulk of the system, decaying as a power of the 
distance I from the surface with: 



mal perturbation in ( 1 . 1 ) transforms as: 



[AK(l)}> 



(1.2) 



AK{i) = 



A 

F 



(i.i) 



Such extended surface perturbations may change the surface 
critical behaviour for an arbitrarily small value of the perturba- 
tion amplitude A. This can be explained [Q] by noticing that, 
under a length scale transformation I' = l/b, the extended ther- 
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V K l h 

so that, comparing both sides of the last equation, the perturba- 
tion amplitude scales as: 

A' = b- K+1/v A. (1.3) 

When the decay is strong enough (k > 1/v), the perturbation 
is irrelevant: it does not affect the surface critical behaviour. 
When k < 1/v, the perturbation is relevant: the decay is suf- 
ficiently slow to modify the surface critical behaviour. In the 
marginal situation, k=1/v, one expects a nonuniversal surface 
critical behaviour. 

Analytical results obtained for the two-dimensional (2d) 
Ising model, with v=l, are in complete agreement with these 
predictions 

i) For k > 1, the scaling dimensions of the surface magneti- 
zation and the surface energy density are x s m = 1/2 and xl = 2, 
respectively, like in the homogeneous semi-infinite system. 

ii) In the marginal case, k = 1, when A is smaller than 
a critical value A c , the surface transition is second-order with 
continuously varying surface scaling dimensions x s m (A) and 
x%(A). For A > A c , the transition is first-order: there is a non- 
vanishing spontaneous surface magnetization at the bulk crit- 
ical point which vanishes above T c since the surface is one- 
dimensional. 

iii) For n < 1, the surface transition is first-order for A > 
and continuous for A < where the surface magnetization 
displays an essential singularity. 



2 



L. Turban: Aperiodic extended surface perturbations in the Ising model 



In the present work we study an aperiodic version of the 
HvL model where the amplitude of the decay is modulated 
according to some aperiodic sequence. More specifically, we 
consider a layered semi-infinite system with constant reduced 
interactions K\ — 0J\ between nearest neighbours in the direc- 
tions parallel to the free surface while the reduced interactions 
in the perpendicular direction are modulated and reads: 



.4, 



(-l) f 'A, 



(1.4) 



where /; is aperiodic and takes the values or 1. The aperi- 
odic sequence is assumed to have equal densities of and 1, 
such that the mean amplitude vanishes. Otherwise, the critical 
behaviour would be the same as for the HvL model with A 
replaced by the average of the modulation amplitude Ai . 

A relevance-irrelevance criterion will be proposed, which 
is a combination of the arguments of Burkhardt and Cordery 
for HvL perturbations [Q| and those of Luck for aperiodic sys- 
tems [0]. 

In Section^, we recall the main results about the fluctuation 
properties of aperiodic sequences generated via substitutions 
and derive a relevance-irrelevance criterion for aperiodic HvL 
perturbations, which is valid to linear order in the perturbation 
amplitude. In Section f|, we examine the critical behaviour of 
the surface magnetization in the aperiodic HvL Ising quantum 
chain. The results of Section ||] are confronted to numerical data 
obtained via finite-size scaling (FSS) for different aperiodic 
sequences in Section f|. The influence of higher-order terms 
on the relevance of the perturbation, the scaling behaviours of 
the first gap and the surface energy are discussed in Section |5|. 
Technical details are given in Appendices A and B. 



2 Fluctuation of the perturbation amplitude 
and relevance-irrelevance criterion 



We consider the aperiodic HvL model defined in (1.4) for a d- 
dimensional semi-infinite system with a bulk correlation length 
exponent v. The modulation of the perturbation amplitude fol- 
lows some aperiodic sequence which can be generated through 
substitution rules [8], 

Working with a finite alphabet A, B, C, ... and starting for 
example with A, an infinite sequence of letters is obtained by 
first substituting a finite word S(A) for A and then iterating 
the process (with B -> S(B), C -> S{C), ...). The required 
aperiodic sequence of digits /; is finally obtained by replacing 
each of the letters by a corresponding finite sequence of digits 
Oandl. 

The information about the properties of the sequence is 
contained into the substition matrix with entries riij giving the 
numbers of letters i in S(j) (i, j = A, B, C, ...). Let V a be the 
right eigenvectors of the substitution matrix and A a the corre- 
sponding eigenvalues. The leading eigenvalue, A\, is related to 
the dilatation factor b of the self-similar sequence and gives its 
length after n substitutions through L n ~ A™. 

The asymptotic densities of letters i = A, B, C, ... are re- 
lated to the components of the eigenvector corresponding to the 



leading eigenvalue through: 



(i) 



Vi(i) 



(2.1) 



which allows a calculation of p^, the asymptotic density of 1 
in the sequence of digits /; . 

The sum of the digits n; — , fa is such that: 

ni - Ipoo ~ l u , (2.2) 

with an amplitude which is log-periodic, i.e., a periodic func- 
tion of ln I j ln b. The exponent 



lnA 1 



(2.3) 



is the wandering exponent of the sequence. When u> < 0(> 0) 
the sequen ce h as bounded (unbounded) fluctuations. 
Using (2.2) and the identity 



■1) 



fk 



1 - 2f k , 



(2.4) 



the amplitude of the aperiodic HvL perturbation in (1.4 ), ave- 
raged at a length scale I, can be written as: 



.4 



k=l 



when the asymptotic density of 1 is poo = 1 /2. 

To analyse the rele vanc e of the perturbation, one may re- 
place Ai by y4/ w_1 in (L4), which leads to the following be- 
haviour under rescaling: 



[AK 2 (l)Y = A \, = b^^ u — — — — 



or, for the perturbation amplitude: 

A' = fj-K+U-l+l/uA 



(2.6) 



(2.7) 



Thus, to the first order in the amplitude, the perturbation is 
relevant (irrelevant) when k < (>)co — 1 + Nonuniversal 
surface critical behaviour is expected in the marginal situation 
where k — uj — 1 + 1/ v. 

For the Ising model in 2d, v=l, and the relevance of the 
aperiodic HvL perturbation is simply governed by the sign of 

LO — K. 

One may notice that Luck's criterion for aperiodic pertur- 
bations [[7]] is recovered when k = 0, which corresponds to an 
aperiodic modulation with a constant amplitude. 



3 Surface magnetization of the Ising 
quantum chain: analytical approach 

In this section, we study the critical behaviour of the surface 
magnetization for the 2c? Ising mod el w ith an aperiodic ex- 
tended surface perturbation given by (1.4). 
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3.1 Ising quantum chain in a transverse field 

We work in the extreme anisotropic limit |J where K\ — ► oo, 
K 2 (l) — > 0, A — > 0, while the ratios 



A = K 2 /Kl 



A/K 2 , 



(3.1) 



are kept fixed. The dual interaction K* — — l/21n(tanh Ki) 
enters into the expression of the unperturbed coupling A. 

In the extreme anisotropic limit, the row-to-row transfer 
operator, T = exp(— 2K*TC), involves the Hamiltonian of the 
quantum Ising chain in a transverse field [P,|lO|]: 



oo 



(3.2) 



i=i 



The as are Pauli spin operators and the inhomogeneous cou- 
plings A, take the HvL form: 



A; 



ft, 



(-i)V 



(3.3) 



where a, is the aperiodic a mpli tude and a is related to original 
parameters as indicated in (3.1 ). A plays the role of an inverse 



temperature and the deviation from the critical point, A c = 1, is 
measured by 



t = 1 - A" 



(3.4) 



so that < t < 1 in the ordered phase. 

Using the Jordan- Wigner transformation [|llj and a cano- 
nical transformation of the fermion operators, the Hamiltonian 
can be put into diagonal form [HOlOj: 



w = E e « 



Vila - g 



where ?yj (?y Q ) creates (destroys) a fermionic excitation. The 
excitation energies e a are obtained by solving the following 
eigenvalue problem: 

Aj-i^-l) + (1 + A?_ x )0 a (O + M a {l + l) = e 2 a (t> a {l) . 

(3.6) 

The physical properties of the system can be expressed in terms 
of the excitations e a and the eigenvectors <p a which are as- 
sumed to be normalized in the following. 



3.2 Surface magnetization 

The surface magnetization m s is given by the matrix element 
(0| erf |cr) , where |0) is the vacuum state and \o) is the low- 
est one-particle state rj\\0). Expressing af in terms of diago- 
nal fermion operators, one obtains m s = 0i(l), which can be 
simply evaluated noticing that the first excitation t\ = in the 
ordered phase [|3j]. This leads to: 



oo J 



S = l 



3 = 1 1=1 



(3.7) 



Using (3.3), the sum can be rewritten as: 

oo 3 



s = i + Y.^ 2]p 3 

3 = 1 



n 

i=i 



1 + a- 



-1)/' 



(3.8) 



Since the critical behaviour is governed by the long distance 
properties of the sytem, when n > 0, one may write: 



lnP J= ^ln 

i=i 

3 



l + a- 



-l)f> 



-iy> 

1F~ 



^- + (-1)^0(1-^) 



•(3.9) 



To go further we need some approximate expression for 
Eqs. ( |2.2| ) and (2.5) suggest that, inside a sum, /, may be re- 
placed by 



/, ~ Poo + C(l) l U 



(3.10) 



where c(l) is a log-periodic amplitude. Making use of the iden- 
tity (2.4), with poo = 1/2, the first sum in the last equation 
of (3.9) can be rewritten as: 



£<-i>* 



i=i 



-2 ]T c(i) r-"- 1 = -2c K (j) £ i^- 1 



1=1 



1=1 



In the following, we assume that 



(3.5) c *ti) 



EUi<i)t 



5^ =1 1 



(3.11) 



(3.12) 



can be replaced by some constant effective value c when j ^> 1 
(see Appendix A for a discussion of the asymptotic properties 
<>!>,. ;./;). 

Finally, when = 1/2 and k > 0, one obtains: 



3 V 



InR 



E 

i=i 



2ca a 2 

+ 



li-w+ K 21 2k 



0{l 



-Zk+uj-1^ 



, (3.13) 



where the size of the leading omitted term follows from a com- 
parison to the first one. Otherwise, when a, ^ 0, an additional 
term of order l~ K gives the leading contribution to the sum 
since u> < 1. As mentioned previously, the critical behaviour 
would then be the same as for the HvL model. 



3.2.1 Relevant perturbations 

We first consider the case of relevant perturbations which, ac- 
cording to the cr iterion of Section ||, corresponds to k<lu. 

In Eq. (3.13) the first term dominates when 1 — u> + n < 2k, 
i.e., when 1 — u < k < ui, which can be satisfied only when 
u! > 1/2. Details about the calculation of the t-dependence of 
m s and its finite-size behaviour at A c are given in Appendix B. 



Fig. 1. Relevance-irrelevance of aperiodic extended surface perturba- 
tions in the Ising quantum chain. The perturbation is relevant in the 
shaded region. To first order in the amplitude a, the perturbation is 
relevant only below the first diagonal. Below the second diagonal and 
the line k = 1/2 the second-order term is relevant and dominates the 
first-order one. On the marginal border line, the magnetization expo- 
nent is a function of a above u = 1/2, of a 2 below, and involves both 
terms at u = 1 /2. The values of uo for the aperiodic sequences used in 
the numerical FSS study are indicated. 



For k = lj > 1/2, the first term of the sum in (3.13) is the 
dominant one. When ca > — 1/4, as shown in Appendix B, this 
term leads to a second-order surface transition where m s van- 
ishes as t^" or, at the critical point, as L~ x ™. Since v = 1, here 
and in what follows, both exponents have the same expression: 

(3 s = x s m = ~(l + 4ca), (3.20) 

When ca < —1/4, the transition is first-order and the approach 
to the critical surface magnetization is governed by the expo- 
nents (3' s — x' s m with 

# = -2&(a) =-l-45a, (3.21) 

where (3 S (a) is the continuation of f3 s in the fi rst-order region. 



For k = (J = 1/2, the two terms in (3.13) are of the same 
order l^ 1 . Thus one obtains a second-order surface transition 
with 

/3 S = i(l + 4m + a 2 ) (3.22) 

for values of a leading to a non-negative expression, i.e., when 
c 2 < 1/4 for any a and when c 2 > 1/4 for a < —2c — A or 
a>-2c+A where Z\=(4c 2 -l) 1/2 . 

When c 2 > 1/4, there is a first-order surface transition on 
the interval —2c— A < a < —2c+A where m s approaches its 
critical value as a power of t or L _1 with an exponent 



P' s =-2j3 s (a) = -l-4ca-a 2 . (3.23) 



L. Turban: Aperiodic extended surface perturbations in the Ising model 



5 



As mentioned above, the second term in (3.13) is the dom 



inant one and leads to a stretched exponential surface magne- 
tization in the region « < 1 — oj, n < 1/2. The marginal line 
k = oj of the linear relevance-irrelevance criterion is actually 
moved to k = 1/2 when oj < 1/2. On this line, the surface 
transition is second-order with (see Appendix B): 



ft 



(1 + a 2 ) 



(3.24) 



The behaviour in the (k, w)-plane is recapitulated in Fig. Ill 
Notice that on the line n = the bulk is homogeneously aperi- 
odic and, according to the Luck criterion in Eq. (2.7), the per- 



turbation is marginal at oj = 0. This marginal behaviour has 
been indeed confirmed by various exact results for aperiodic 
Ising models 00). 

In Fig. [j] this point belongs to the domain of relevant pertur- 
bations due to terms of order a 2 . This only means that the non- 
decaying aperiodicity changes the critical coupling from A c = 1 
to A c = (1 — q 2 )~ 1 / 2 [^6|]. At the unperturbed fixed point value 
of A, the aperiodic system is in its disordered phase. Thus due 
to terms of order a 2 , the aperiodicity induces a flow towards 
the trivial fixed point where the coupling vanishes. 

One may notice that the methods used above do not apply to 
the case of ho mmog eneo us ape riodic systems: the correspond- 
ing equations, (3.18) and (3.19), are valid only when k > 0. 



4 Surface magnetization of the Ising 
quantum chain: finite-size-scaling study 

The results obtained in the last section rely on the v alidity of 
the assumptions used to transform Eq. (3.9) into Eq. (3.13). In 



this Section these results are checked through a numerical FSS 
study of the critical surface magnetization for different systems 
corresponding to the different regions of the (k, w)-plane in 
Fig. pi We also evaluate numerically the effective amplitude c 
which is needed when the term linear in a contributes to the 
critical behaviour. 



4.1 Aperiodic sequences 

We used aperiodic sequences with the same asymptotic density 
Poo = 1/2 leading to different values of the wandering exponent 
oj and also to different behaviours for the log-periodic ampli- 
tude c K (j). 

The Rudin-Shapiro (RS) sequence [|l7[ follows from sub- 
stitutions on the letters A, B, C, D, with S(A) = AB, S(B) = 
AC, 5(C) = DB, S(D) = DC. The substitution matrix has 
eigenvalues Ai — 2 and A% = y2 so that oj = 1/2. The different 
letters have the same asymptotic density poo(t) = 1/4 (i = A, 
B, C, D). Each letter in the sequence corresponds to a pair 
of digits A = 00, B = 01, C = 10 and D = 11, which gives 
Poo = 1/2. 

The same substitutions on the letters are used to generate 
the RS8 sequence which is obtained by replacing and 1 by 
0000 and 1111 in the the RS sequence. The substitution matrix 
being the same, both the wandering exponent and the asymp- 
totic density keep their RS values. The new correspondance 



6.0 



4.0 



2.0 



0.0 



-3.0 




a 



Fig. 2. Amplitude C(a) for the TM1 sequence with k — 1/2 deduced 
from FSS at criticality with L n = G n , n = 1, 7. The surface transition 
is continuous for a < and first-order for a > 0. 




-0.1 



-0.2 - 



-0.3 



-0.4 



4.0 6.0 8.0 

Inj 



Fig. 3. Amplitude c K {j) as a function of lnj for the TM1 sequence 
with k = 1/2. Asymptotically c K (j) oscillates log-periodically be- 
tween the two values indicated on the left axis. 



between letters and digits only affect the behaviour of the log- 
periodic amplitude c K (j). 

The RS and RS8 sequences are self-similar under dilatation 
by a factor b = 4. 

Other values of the wandering exponent have been obtained 
through decoration of the Thue-Morse sequence which follows 
from the substitutions 5(0) = 01 and 5(1) = 10 [@. 

The TM1 sequence is generated through 5(0) = 010000 
and 5(1) = 101111. The two eigenvalues of the substitution 
matrix are A x = 6 and A 2 = 4 so that oj = In 4/ In 6 ~ 0.7737. 
Due to the invariance under the exchange of and 1 , the two 
digits have the same asymptotic density = 1 /2. 

Finally 5(0) = 010100 and 5(1) = 101011 lead to the TM2 
sequence. Its substitution matrix has eigenvalues A\ = 6 and 
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0.5 




20 
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10 



0.0 1 — ■ — 1 — ^ — — 1 — 1 — ■ — 1 — ■ — 1 
-1.0 -0.5 0.0 0.5 1.0 1.5 2.0 

a 

Fig. 4. Amplitude C(a) for the TM 1 sequence with k, — 1 — w deduced 
from FSS at criticality with L n = 6™, n = 1, 5. The surface transition 
is first-order in the intermediate region where C(a) vanishes and con- 
tinuous elsewhere. 



o 1 — ■ — 1 — ■ — 1 — — 1 — ■ — 1 — ■ — 1 

-3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 
a 

Fig. 6. Amplitude C'(a) for the TM2 sequence with k — u deduced 
from FSS at criticality with L n — 6", n= 1, 7. The surface magnetiza- 
tion vanishes continuously with a parabolic amplitude in the stretched 
exponential. 
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1.0 



0.5 



0.0 



-1.2 -0.8 -0.4 



0.0 

a 




Fig. 5. Amplitude C(a) for the RS sequence with k = 0.3 deduced 
from FSS at criticality with L n =4 n ,n — 1,9. The surface magnetiza- 
tion vanishes continuously with a parabolic amplitude in the stretched 
exponential. The solid (dashed) line corresponds to a correction-to- 
scaling exponent 0.5 (0.01) in the BST extrapolation. 

A 2 = 2, which gives ui — In2/ln6 ~ 0.3869. For the same 
reason as above = 1/2. 

The TM1 and TM2 sequences are self-similar under dilata- 
tion by 6 = 6. 



4.2 Relevant perturbations 



In the case of relevant perturbations, according to ( B.4 ), ln m s>c 
is linear in L T with a slope —C(a) when the transition is con- 
tinuous. We studied the size-dependence of lnm s c for sizes of 
the form L n ~ b n . Approximants C n (a) for the amplitude, de- 
duced from two-point fits for successive sizes L n , L n+ i, were 



extrapolated using the BST algorithm [|l8|]. When the transition 
is first-order, the leading contribution is independent of the size 
of the system and C(a) vanishes. 

With the TM1 sequence at k= 1/2, the critical behaviour 
is governed by a term of order a (see Fig. |l]). The extrapolated 
amplitude C(a) is shown in Fig. |[ It displays the expected 
linear variation with a in the region a < where the transition 
is continuous. The slope leads to the effective amplitude c given 
in Table 1 . 

Since k < uj the amplitude c K (j), which is shown in Fig. pi 
remains log-periodic at infinity (see Appendix A). For a < 
the main contribution to PJ 2 in ( jB.6| ) comes, at large j-values, 
from the absolute minima of c K (j) and c in Table 1 corresponds 
to the extrapolated value (c K ) m i n . 

For the TM1 sequence with k = 1 — uj the surface magne- 
tization at criticality is given by ( 3.17 ). In Fig. [j], this system 
is on the solid line inside the domain of relevant perturbations 
where both a and a 2 contribute to the critical behaviour. 

The FSS study leads to the extrapolated amplitude C(a) 
shown in Fig. [i| We were limited to sizes up to 6 5 because the 
exponent of L in the stretched exponential is two times larger 
than before and m SjC (L) becomes quite small at large size. The 
corresponding parameters are given in Table 1 . 

As above c K (j) remains asymptotically oscillating, here be- 
tween (c K ) min = -0.34(1) and (c K ) max = -0.055(1). 

When a < the transition is continuous and C (a) displays 
the expected parabolic variation. The coefficient of the linear 
term leads to an effective amplitude c in good agreement with 
the value of (c K ) m i n which gives the main contribution to P^ 2 
in ( |B~7l ). 

When a > the transition is first-order below a critical 
value a c . In the region a > a c where the transition is conti- 
nuous, c is close to (c K ) max as expected but the agreement is 
poor for the coefficient of the quadratic term. In both cases the 
sizes used in the FSS study are too small to obtain truly reliable 
estimates. 
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sequence 


TM1 


TM1 (o<0) 


TM1 (a>a c ) 


RS 


TM2 


UJ 


0.7737 


0.7737 


0.7737 


1/2 


0.3869 


K 


1/2 


l-U) 


1 - u 


0.3 




c from c K (j) 


-0.272(7) 


-0.34(1) 


-0.055(1) 






c from C(a) 


-0.271(5) 


-0.332(3) 


-0.04(3) 






72 (expected) 




0.9134 


0.9134 


5/4 


2.20951 


72 (numerical) 




0.92(2) 


0.6(1) 


1.251(1) 


2.209(1) 



Table 1. Parameters of the amplitude C(a) in 
the stretched exponential deduced from the FSS 
study of the surface magnetization at criticali- 
ty. 72 is the coefficient of the quadratic term in 

C(a). 



2.5 



2.0 
1.5 
1.0 
0.5 




0.0 1 — ■ — 1 1 — 1 — ■ — — 1 — ■ — 1 

-5.0 -3.0 -1.0 1.0 3.0 5.0 
a 

Fig. 7. Surface magnetization exponents for the TM1 sequence with 
k — u). The solid line corresponds to (3 B at the second-order surface 
transition and the dashed line, in the first-order region on the right, to 
f3g which controls the approach to the critical magnetization. 



-0.1 




-0.4 u ■ 1 ■ 1 ■ 1 ■ 1 

0.0 2.0 4.0 6.0 8.0 

Inj 



Fig. 8. Amplitude c K {j) as a function of lnj for the TM1 sequence 
with k = u). Asymptotically c K (j) converges towards the effective 
amplitude c indicated on the left axis. 



The RS sequence with k = 0.3 leads to a relevant pertur- 
bation for which th e term in a 2 governs the critical behaviour 
as indicated in Eq. (3.19). The extrapolated amplitude C(a) is 
shown in Fig. || and the parameters given in Table 1. 

For the solid line, obtained with a correction-to-scaling ex- 
ponent equal to 0.5 in the BST extrapolation, the parabola is 
not centered on a — 0: there is a weak linear contribution to 
the amplitude. With a correction-to-scaling exponent equal to 
0.01 (dashed line in Fig. ||) the extrapolation is less stable but 
the coefficient of the linear term is reduced from —0.15 to 
—0.09. Thus we suspect that the unexpected linear contribu- 
tion to C(a) is a correction-to-scaling effect. One may notice 
that a power of L in front of the stretched exponential leads to 
a logarithmic correction. 

Next we consider the TM2 sequence with k = oj. The cor- 
responding point in the (n, cj)-plane belongs to the dashed line 
in Fig. [j]. It leads to a perturbation which is marginal to linear 
order in a but becomes relevant due to the term of order a 2 . 
The extrapolated amplitude is shown in Fig. [| 

The amplitude C(a) is still given by Eq. ( [3.19 ) and dis- 
plays a parabolic behaviour. The coefficient 72 in Table 1 is in 
goog agreement with the expected value. The same parabolic 
variation was obtained for the TM2 sequence with n = 0.3 and 
K=QA. 



CO. 



4.0 



3.0 



2.0 - 



1.0 



0.0 1 1 1 1 1 1 1 1 1 1 — 

-2.0 -1.0 0.0 1.0 2.0 3.0 

a 

Fig. 9. Surface magnetization exponent f3 a for the RS sequence with 
h — uj. The surface transition is always second-order. 



4.3 Marginal perturbations 

The critical surface magnetization has been calculated using 
the same chain sizes as in the relevant situation (up to n = 7 
for TM1). Estimates for the exponent f3 s are obtained via two- 
point fits of ln[m SiC (£n)] versus lnL„. The two-point approx- 
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sequence 



TM1 



RS 



RS8 



TM2 



LO 


0.7737 


1/2 


1/2 


0.3869 


K 


a; 


a; 


LO 


1/2 


cfrom c K (j) 


-0.199(1) 


-0.377(4) 


-0.75(1) 




c from /3 S (a) 
c from /3s (a) 


-0.1998(1) 


-0.3763(1) 


-0.7526(2) 




-0.19(1) 




-0.7525(15) 




/3o 


0.49998(2) 


0.4997(3) 


0.4999(5) 


0.502(5) 


& 




0.50007(7) 


0.5000(1) 


0.499(1) 




-0.99(2) 




-0.999(1) 










-0.999(1) 





Table 2. Parameters of the magnetic exponents /3 s (a) and 
f3' s (a) deduced from the FSS study of the surface magneti- 
zation at criticality. The expected value for the constant term 
/3o (/3q) and the coefficient of the quadratic term /32 in 
/3 s (a)isl/2(-l). 



-0.3 




Fig. 10. Amplitude c K (j) as a function of lnj for the RS sequence 
with k = lo. Asymptotically c K (j) converges towards the effective 
amplitude c indicated on the left axis. 



imants are extrapolated using the BST algorithm. In the first- 
order regions, the regular contribution associated with the crit- 
ical magnetization m s c (oo) is eliminated using the differences 



Am s ^{L r , 



z(L n +i) and the usual procedure 



is applied to calculate the exponent f3' s , although with one size 
less. 

As a first example of marginal behaviour we consider the 
TM1 sequence with k = lo. The extrapolated exponent values 
are shown in Fig. ^| 



In agreement with Eqs. ( 3.2C ) and ( 3.21 ), both exponents 
vary linearly with a. Due to the finite sizes used, the singu- 
larities remain rounded near a c where the surface transition 
changes from second to first order. 

As shown in Fig. |[ c K (j) here converges towards the ef- 
fective amplitude c. It is compared to the values deduced from 
the slopes of the surface exponents in Table 2. The precision on 
the parameters deduced from f3' s (a) is lower since the number 
of points is reduced by one in the extrapolation. 

The RS sequence with n = lo = 1/2 leads to a marginal 
behaviour with linear and quadratic contributions to the expo- 
nents. The exponent (3 S is shown in Fig. ^. The corresponding 
parameters are given in T able 2. The variation is parabolic in 
agreement with Eq. ( 3.22 ) and the absence of a first-order re- 
gion is linked to the small value of |c|. The log-periodic am- 



5.0 ry 

4.0 - 

3.0 - 

2.0 - 
1.0 
0.0 



-2.0 -1.0 0.0 



1.0 2.0 
a 



3.0 4.0 5.0 



Fig. 11. Surface magnetization exponents for the RS8 sequence with 
n = lo. The solid line corresponds to j3 a at the second-order surface 
transition and the dashed line to /3' s which controls the approach to the 
critical magnetization. The chain sizes used in the FSS study are of 
the form L n = 2 x 4 n , n= 1, 8, with one size less for f3' B . 



plitude c K (j) shown in Fig. [10| converges to c = —0.377(4) 
whereas the limiting value for the occurence of a first-order 
transition is \c\ = 1/2. 

With the RS8 sequence at k = lo = 1/2 we obtain once 
more a marginal perturbation. But now c K (j) converges to a 
value allowing for the occurence of a first-order transition. The 
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exponents (3 S and j3' s are shown in Fig. 

The surface transition is first-order in the central region and 
second-order outside, with a parabolic variation of the expo- 
nents in both cases , in a greem ent with the analytical expres- 
sions in Eqs. ( [3.22 ) and (3.23). The accuracy of the extrapola- 
tion is reduced near the singularities on the borders. The coeffi- 
cients given in Table 2 are in good agreement with the expected 
ones. 

Finally we have studied the TM2 sequence with n= 1/2 as 
an example of a system for which the marginal behaviour is in- 
duced by the term of order a 2 . In this case the surface transition 
is always second order. The exponent B s in Fig. |l2| displays the 
parabolic behaviour obtained in Eq. (3.24). Here too the coef- 



ficients, given in Table 2, are close to the expected values. 



Fig. 13. Relevance-irrelevance of aperiodic extended surface pertur- 
bations in a system with bulk correlation length exponent v. The per- 
turbation is relevant in the shaded region. The order of the dominant 
perturbation terms is indicated, assuming that off-diagonal terms do 
not contribute. 



One may notice that such cross-terms do not enter into the cal- 
culation of the surface magnetization of the 2d Ising model. 

The diagonal contributions to involve pairs of sites 
belonging to the same layer and read: 

-^ 2) = ^EE^Mni) ; (5-4) 



where is the surface of the layers. The density under the 
sums has dimension 2d — 1 whereas the correlation function 
has dimension 2x c , where x c is the scaling dimension of the 
bulk energy density, thus one obtains: 



( A ' 2 )d _ h 2d-l-2x e ( A2 )d _ ,-1+2/u i A2 )d 

1' 2k 1 2k 1 2k 



(5.5) 



and the diagonal second-order amplitude transforms as: 

(A' 2 ) d =b- 1 - 2K+2/ »{A 2 ) d . (5.6) 

It is a relevant (irrelevant) variable when k < (>) — 1/2 + 1/v 
and a marginal one when k = — 1/2 + 1/v. 

The scaling dimension U(A 2 ) d of (A 2 )d has to be compared 



to the dimension jja of A given in (2.7) to see which term 
governs the critical behaviour when both are relevant. When 
U(A 2 ) d > UA, i-e-, when k < —u) + 1/v, the second-order term 
dominates and one expects an A 2 -dependence of the ampli- 
tudes in the stretched exponentials. On the line k = —oj + 1/v 
when k<— 1/2 + 1/v, linear and quadratic terms contribute 
together. 

A summary of the relevance-irrelevance in the (k, w)-plane 
is given in Fig. [l3|. It is in agreement with our findings of Sec- 
tions |^ and ^ for the surface magnetization of the 2c? Ising 



10 



L. Turban: Aperiodic extended surface perturbations in the Ising model 



model. The order of the dominant contributions could be mo- 
dified by the discarded off-diagonal second- or higher-order 
terms. As mentioned earlier for the 2d Ising model, there is 
a shift of the bulk critical coupling for the homogeneous aperi- 
odic system on the line k = 0. Along this line the perturbation 
is irrelevant below u = 1 — 1/v. 

One may notice that the model shows some kind of univer- 
sal behaviour. Both the amplitude C(a) in the stretched expo- 
nential and the surface magnetic exponent j3 s (a) are indepen- 
dent of the aperiodic sequence, provided the value of w is such 
that the system belongs to the region k < 1/2, n < 1 — u> of 
Fig. [l] where the quadratic contribution dominates. 

Let us now briefly discuss the behaviour of the first gap and 
the surface energy. 

On ajinite critical Ising chain, the first gap is known to 
scale as 



(5.7) 



where m s c is the critical magnetization on the second surface 



and Pl is the product of the couplings defined in Eq. (3.8) 



When k > the bulk is unperturbed and, asymptotically, the 
second surface displays an ordinary surface transition. Thus the 
scaling dimension of m s is x^ n = 1/2. 

For relevant perturbations, with the notations of Appendix 
B, we have: 



Pr. 



(5.8) 



Wh en th e surface transition is continuous m SiC (L), given 



in Eq. (B.4), behaves as Pl- It follows that the first gap van- 
ishes as a power of L. One expects in this case the unperturbed 
behaviour £i(L)~L _1 . 

When the surface transition is first-order, the leading term 
in m StC (L) is a constant and there is no more compensation. 
The first gap is anomalous, it vanishes with an essential singu- 
larity. This behaviour is linked to the localization of the cor- 
responding eigenvector cf> 1 , which itself is reponsible for the 
finite weight on the first component <f>i(l), leading to a non- 
vanishing surface magnetization [p3|. 

The scaling dimension a;| of the surface energy can be de- 
duced from the finite-size behaviour e s = (0|crf |e) where the 
state |e) = ^jf^lO) is the lowest two-particle excitated state. 
This matrix element can be written as [ttgh: 



e s = (e 2 -ei)<Ml)0 2 (l). 



(5.9) 



For relevant perturbations one expects an essential singularity 
since at least ^(l) should behave in this way. 

In the case of a marginal perturbation, with the notations of 
Appendix B, one may write: 



P; 



L B/2 _ £-l/2+,3 s (a) 



(5.10) 



where (3 B (a) is the surface magnetic exponent f3 s when the 
transition is second-order or its continuation to negative values 
when the transition is first-order, i.e., when /3 S = 0. 

When the transition is second-order Eqs. ( [5.7| ) and ( 5.1C ) 
lead to: 



(5.11) 



i.e., to the unperturbed behaviour as above for relevant pertur- 
bations when the transition is continuous. 

For the surface energy, Eq. ( ^^ ) leads to the scaling dimen- 
sion 



= 1 + 2/3 s 



(5.12) 



if one assumes that 02(1) scale as L _/3s too. This is known 
to be true either for the marginal HvL model [ pp] ] or for the 
aperiodic version of the same model [ |ll^ ] where the couplings 
are modulated according to the Fredholm sequence Jl7|]. 

When the transition is first-order, due to the localization of 
4> x , the scaling of the first excitation is anomalous [plj]: 



ei(L) ~ L- 1+ *(°> 



(5.13) 



It decays faster than higher excitations since s (a) < 0. 

For the surface energy, we conjecture the following beha- 
viour: 



i- A(o). 



(5.14) 



The factor containing the excitations in Eq. ( |5.9| ) is dominated 
by €2 which vanishes as L^ 1 . Furthermore we assumed that, 

like in Refs. |jl|,|§, cfc(l) scales as L~ P ' J2 - 

As for the HvL model, in the regime of first-order transi- 
tion, the anomalous scaling of the first gap leads to an exponent 
asymmetry For example, in the disordered phase, the expo- 
nent of the correlation length ^ , along the surface of the semi- 
infinite system, is governed by the first gap and v\\=l — f3 s (a). 
In the ordered phase the first excitation vanishes and £|| ~ e^ 1 
so that j/jj — 1, like in the unperturbed system. 

To conclude, let us mention that the case of random surface 
extended perturbations, which has fluctuation properties simi- 
lar to the aperiodic case with u> = 1/2 J22|], is currently under 
study. 



I thank Ferenc Igloi and Dragi Karevski for constructive comments 
and a long collaboration on aperiodic systems. 



Appendix A: Log-periodic functions 



The log -periodic function c(l) defined in Eq. ( 3.10 ) is such that 
c(b n l) = c(l) where b is the discrete dilatation factor of the 
aperiodic sequence. It can be generally written as a Fourier ex- 
pansion, 



c(l) = c + ^2 Ck cos[0fc(Z)] , 



fc=i 



mo 



so that, in (3.12): 



cos 8 k (I) 



k=l 



-K-l 



(A.l) 



(A.2) 
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We are interested in the behaviour of c K (j) at large j. 
• When k<u>, let us write: 

oo 

c K (j) = c + ^2 c k $t[I k (j)] , 



(A.3) 



fe=i 



where: 



h(j) 



~ e 



EO /w-K-l+i27rfc/lnb 
, 1=1 1 

ci j; ;a;-«-l+i27rfc/ln6 
i<Pk ±0 



(A.4) 



so that: 



cos[6» fc (j)] + a fc sin[0 fc (j)] 



1 + a 2 



2irk 



Oik = 



(lj — k) In b 



(A.5) 



It follows that c K (j) oscillates log-periodically around cq. The 
Fourier coefficients Ck are divided by k whe n k is large or when 
k is close to u). The effective constant c in ( 3.13 ) can be taken 
as the value of c K (j) which gives the main contribution to the 
function in which it enters. 



• When k — uj, replacing the sums over I in ( \.2) by inte- 
grals, the change of variable u = ln I leads to: 



°° c k f lnj ( u 

c K (j) ~ c + — / du cos [2nk— - 



fk 



(A.6) 



Thus liirij—xx, c K (j) = Co and c is the constant term in the 
Fourier expansion of c(/). 



• When k > lo, according to ( A.4 ), the asymptotic expres- 
sion of c K (j) can be written in terms of £ functions with: 



ifc(oo) 



jip k C(l + « - w - i27r/c/ In 6) 
C(1 + k- w) 



(A.7) 



Thus, in this case too, c K (j) tends to a well-defined limiting 
value giving c. 



Appendix B: Amplitudes and exponents 

We study successively the temperature-dependence of the sur- 
face magnetization near the critical point as well as its size- 
dependence at criticality. 

We first consider the case of relevant perturbations where 
the leading contribution to I n Pj is some positive power of j. 
The sum S = m~ 2 in Eq. (3.8) can be then replaced by an 
integral of the form [H[23J 



p OO pOC 

S(t) ~ / dj \- 2 Jpr 2 ~ / dj exp 
Jo Jo 



(-tj + Bf) , 



(B.l) 



where we used the definition of t given in (3.4). For the finite- 
size behaviour at the critical point, t = 0, the sum is cut off at L 
so that: 



S C (L) 



[ djPT 2 - f djexp(Bf), (B.2) 
J o Jo 



When B > 0, the integral can be evaluated using Laplace's 
method when < r < 1 and, up to a power law prefactor, one 
obtains: 



m s (t) ~ exp 



1-T 

2r 



V(1-t) +-r/(l-r) 



(B.3) 



The main contribution to S C (L) comes from the vicinity of the 
upper limit. Expanding the argument of the exponential around 
L leads to: 



s ,c(£) ~ exp 



— —B L T 
2 



(B.4) 



When B < 0, expanding cxp(— ij) in (B.l) and integrating 



term by term 


gives: 




m s (i) ~ 


"tISI 1 /-" 


1/2 , 







r(2/r) 



2r(l/r)| J B| 1 A 



(B.5) 



Let us now look for the values of B and t when k varies. 
• W hen 1 — lo < k < uo, the term which is linear in a domi- 
nates in (3.13), so that: 



P: 



exp 



4ca 



-r 



(B.6) 



which leads to the expressions given in E qs. (B.14 ) and (3.15). 

• When 1 — lo = n < u>, both terms in ( 3 . 1 3 ) are of the same 
order and 



Pr 2 ~ exp 



4ca + a 2 . x _ 2k 
1 - 2k 



-J 



(B.7) 



leading to ( |3.16[ ) and ( p.l7| ). 

• When k < m in(l — u, 1/2), the term in a 2 governs the 



behaviour of ( 3 . 1 3 ), so that 



p; 



exp 



1 - 2k 



1— 2k 



(B.8) 



from which Eqs . (3.18) and (3.19) follow. 

Next we consider the case of marginal perturb ation s where 
P 3 behaves as?'" 5 / 2 . When B > -1, S in Eq. (Q can be 
rewritten as Jl3|] 



/•OO PC 

S(t)c / <1,A 2 T/ - / 
Jo Jo 



or, at the critical point, 



S C (L) 



djP- 



djj B e-V, 



djj 1 



(B.9) 



(B.10) 
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For the ^-dependence, one obtains: 



m s (t) 



-l-B 



du u e 



B„-u 



-1/2 



whereas: 



[riB + i)}- 1 / 2 ^ 1 ^ 2 , 



3 (Z)~(1+B) 1/2 £~ (1+B)/2 



(B.ll) 



(B.12) 



When B< — 1, the integrals in (B.9) and ( ft. 10 ) diverge at their 
lower limits. The main contribution, coming from small values 
of j, must be treated more carefully. For this purpose let us split 
S into two parts as: 



s(t) = i + £ r i*i+5>-«-i)j 



-l-B 



(B.13) 



The first sum gives whereas the second sum can be 

transformed using the Euler-MacLaurin summation formula: 



£ /(?',*) = / d i/0'. f )+ / d J 

Jl 



The change of variable u = tj leads to: 



f(j,t) + 0(t). 

(B.14) 



S(t) = m-f + i 



-2 , ,|B|-1 



du (e~ u - l) U -^ 



(B.15) 



where m s . c = [1+£(|.B|)] _1 ' 2 . The i-dependence of the first in- 
tegral is obtained through an expansion of the exponential and 
the second is a constant. Collecting the different terms gives 



S(t) = m~c + const + 0(t) 



(B.16) 



and 



m s (t) = m SiC + const + 0(t) . (B.17) 

For the size-dependence at criticality, one may write: 

L 



OO 



^ p OO 

i+Er |BL dj r w + o( L -W) 

3=1 L 



+ 0(L-\ B \). 



so that: 



m s , c {L) = m SjC + 



2(151-1) 



L 



-\B\+1 



(B.18) 



(B.19) 



Finally, we identify the expression of the exponent B for 
K-values where a marginal behaviour is obta ined. 

• When k = u) > 1/2, the linear term in ( 3.13 ) is the domi- 
nant one and yields: 



(B.20) 



from which the exponents in (3.2C) and (3.21 ) follow. 

• Wh en k = oj — 1/2, the linear and quadratic terms in 
Eq. (3.13) contribute, so that 



p — 2 ^ -ica+a 
3 ~ ■' 



(B.21) 



leading to ( p.22[ ) and ( [3.231 ). 

• When k = 1/2 and u> < 1/2, the leading contribution 



in (3.13) is the quadratic one, 



(B.22) 



which gives the exponent in Eq. (3.24). 
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